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1. Introduction 

Interest in strongly coupled gauge-fermion models has skyrocketed in recent years. Lattice 
methods offer the most reliable tools, and in the Lattice 2011 Conference there were over two 
dozen contributed talks discussing these systems. Many of the talks considered SU(3) gauge fields 
with 12 flavors of fundamental fermions, one of the most studied yet still controversial systems. 

2-loop perturbation theory predicts conformal behavior with Nf > 8 flavors, while summing 
up rainbow diagrams for a more reliable estimate puts the conformal window just below 12 fla- 
vors. Early numerical studies used the Schrodinger functional approach with unimproved staggered 
fermion lattice action to numerically calculate the renormalization group j3 function and concluded 
that the theory has an infrared fixed point (IRFP)[1, 2]. Ref.[3] considered the system at finite 
temperature and reached similar conclusion. On the other hand spectral quantities appeared to be 
more consistent with a chirally broken system[4, 5, 6]. The first studies using Monte Carlo renor- 
malization group (MCRG) techniques were not able to push deep enough into the strong coupling 
and remained inconclusive[7, 8]. 

Recently a large scale study showed that an extensive data set of spectral quantities prefers 
the chirally broken interpretation [9], but other groups interpreted the same data as compatible 
with conformal behavior[10, 11]. In a recent work I revisited the 12 flavor SU(3) system using 
MCRG methods. Due to two improvements, one in the lattice action, the other in analyzing the 
MCRG data, I was able to cover a wider coupling range and demonstrate that the renormalization 
group /3 function develops a zero, signaling the existence of an infrared fixed point and infrared 
conformality[12]. Other recent works include Refs. [13, 14], reporting conformal infrared be- 
havior based on the running of a renormalized coupling defined through twisted Polyakov loops. 
Two other works studied the phase diagram of the 12 flavor system [15, 16]. Both studies find 
three phases and indication of two bulk phase transitions, suggesting conformal behavior, but these 
works interpret the nature of the phase transitions quite differently. 

Undoubtedly more work is needed to understand why different methods point to different 
conclusions and to find a consistent picture of the Nf = 12 SU(3) system. The MCRG method is 
just one of the possible approaches to study the fixed point structure of the model, but the necessity 
to use an improved gauge action discussed in Ref. [12] is quite general and applicable to all other 
investigations. In the following I summarize the results presented in Ref. [12] and include some 
new results further supporting the conformal scenario of the Nf = 12 model. 

2. Phase diagram of a system of 12 fermions with fundamental-adjoint gauge action: 

The basic observation that makes the modified action necessary is the existence of an ultravi- 
olet fixed point due to strong coupling lattice artifacts. It is well known that the pure gauge SU(2) 
and SU(3) theories have a first order phase transition line in the fundamental and adjoint plaquette 
action space[17, 18]. This line ends in a second order point that corresponds to a (most likely 
trivial) ultraviolet fixed point (UVFP). While the first order phase transition and the new UVFP 
are lattice artifacts and independent of the Gaussian FP and the continuum limit defined there, their 
existence can strongly influence, even completely change, the scaling behavior of the lattice model. 
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For example Ref. [19] studied the scaling of several observables of the SU(3) fundamental- 
adjoint pure gauge system and showed very large scaling violations, even lack of scaling, at cou- 
plings near the extension of the first order phase transition line. As expected, the scaling violations 
decrease with negative adjoint terms in the action, i.e. farther from the second order endpoint. 
Refs. [20, 21] studied the RG flow lines in the pure gauge SU(2) system and found that near the 
extension of the first order phase transition line along the fundamental plaquette action the RG 
flows away from the UVFP-2 and turns around sharply at negative adjoint coupling. This again 
indicates that in this region the system is strongly influenced by the FP associated with the second 
order phase transition. In a recent work [22] we studied the pure gauge fundamental-adjoint SU(2) 
system with the 2-lattice matching MCRG method. Our results show that near the first order line 
and its extension towards negative adjoint couplings the MCRG matching method breaks down, 
the system is no longer in the basin of attraction of the perturbative FR 

When two UVFPs exist numerical simulations have to stay in the vicinity of either one of them 
to describe the corresponding continuum physics. The basin of attraction of the new UVFP is at 
fairly strong coupling, and present day lattice simulations with 2 or 2+1 flavors are sufficiently far 
from it. This, however, might not be the case in many fermion systems where interesting physics 
is expected to occur at strong gauge coupling. 

Large number of fermions could change the phase diagram of the pure gauge system so the 
first step is to study the phase diagram of the fundamental-adjoint plaquette gauge action with 12 
fermions. I used nHYP smeared staggered fermions [23] and measured the plaquette, the specific 
heat through the derivative of the plaquette, the Polyakov line and the chiral susceptibility on 8 3 x 4, 
8 4 , 12 3 x 4 and 12 3 x 6 lattices. The specific heat gave a very clear signal for a first order transition, 
continuing along a crossover line, as indicated by the solid and dashed red lines in Figure 1. This 
phase transition/crossover has no dependence on the temporal lattice size, it is a bulk feature of the 
system. The phase diagram of the 12 flavor system looks very similar to the pure gauge one. 

The horizontal blue line in Figure 1 at jSa = indicates the region I studied in Ref. [8]. I found 
that matching became impossible at stronger couplings, to the left of the blue line, most likely due 
to the nearby crossover region. The leading order perturbative relation between the gauge coupling 
and the lattice couplings is 



This suggests that the coupling (j3f,j3a)=(5. 0,0.0) corresponds, at least perturbatively, to (fip , /3a)=(10.0, 
-2.5). The latter point is quite far from the crossover along the Pa/Pf = —0.25 action line, indi- 
cated by the second blue line in Figure 1 . If the basin of attraction of the perturbative FP is limited 
by the first order/crossover line and Eq. 2.1 is any indication of constant physics lines, than along 
Pa/Pf = —0.25 one can reach considerably stronger couplings than with the Pa =0 fundamental 
action. Finally, the green line in the figure corresponds to Pa/Pf = —0.50, the limit where the 
adjoint plaquette overtakes the fundamental one and flips the system into a new universality class. 

Lattice simulations that want to investigate the properties of the Gaussian FP or the corre- 
sponding IRFP, if it exists, should stay away from the any spurious FP. This is true for any kind of 
investigations, be it Schroedinger functional or other running coupling calculations, spectral studies 
or finite temperature investigations. 




Pf(1+2^). 



(2.1) 



3 



MCRG study of 12 flavors 



Anna Hasenfratz 



6 



GO. 




_4 I I I I I I I I I I I I I I I I I I I I 

2 4 6 8 



Figure 1: The red curve shows the approximate location of the phase transition /crossover in the 
fundamental-adjoint plane for the Nf = 12 flavor system. 
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Figure 2: The optimized Aj3 at j5p = 6.5 with the action Pa/Pf = —0.25 . The red circles are the finite 
volume corrected predictions, the blue diamonds have no finite volume correction in the optimization step. 
The left side of the figure shows results after comparing blocking steps ri/, = 3/2/1 on different volumes, 
while right side is the result after one more blocking step. The data points are labeled by the final blocked 
lattice size Lj,. 
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3. MCRG results with a fundamental-adjoint gauge action: 

The 2-lattice matching MCRG method is a powerful tool to numerically calculate the bare step 
scaling function, the discretized lattice analogue of the RG function [7, 8]. The method relies on 
matching observables after RG blocking steps. Simulations do not have to be performed in volumes 
with lattice size comparable or larger than the correlation length, as most of the finite size effects are 
cancelled by comparing blocked observables measured on identical blocked volumes. The success 
of the method relies on the optimization of a blocking parameter that is set such that it assures the 
FP and its renormalized trajectory of the system is reached in a minimum number of blocking steps. 
It has not been recognized previously that this optimization procedure also introduces finite size 
effects. These effects, while much smaller than the finite size effects coming from the matching 
of observables, can be important in systems where the step scaling function is very small. In this 
section I describe how to minimize this secondary finite volume effects. 

Under repeated renormalization group blocking the action first flows towards the renormalized 
trajectory (RT), than along it, away from the UVFR If the RG flow lines originating at j8 and f5' hit 
the RT at the same points but require nj, and ny — l blocking steps to do so, the conelation lengths 
at /3 and /3' differ by the scale of the block transformation, 2 in the present work. There are two 
steps necessary to achieve matching: 

• Matching: Observables measured on n\, times blocked configurations generated at j8 have to 
match observables measured on — 1 times blocked configurations generated at jS' if the 
blocked actions are identical. Since blocking reduces the lattice size by a factor of two, in 
order to minimize finite size effects it is best to do the simulations on twice as large lattices at 
j8 than at /3' so the final measurements are done on the same volume. The shift in the gauge 
coupling is defined as Aj8^(j8;?i/,,L/,) = j8 — /3' if 



where {G) denotes the expectation value of some short distance operator and Lb is the volume 
after ny and rib — I blocking steps (the same for both sides of Eq. 3.1). 

• Optimization: The quantity A/3 defined in the previous step can differ significantly from the 
step scaling function Sb if the RG flow does not reach the RT in rib — 1 steps. Most RG 
block transformations have a free parameter, usually denoted by 05, that can be optimized 
to minimize the number of RG steps needed to reach the RT. The optimized parameter is 
defined as the one where consecutive blocking steps predict the same shift, 



Previous studies did not take the volume dependence of A/3<? into account and usually satis- 
fied Eq. 3.2 on different volumes. The error introduced this way is small, but is important 
when Afi^ itself is small. Note that matching according to Eq. 3.2 requires simulations on 
three different volumes. 

The procedure can be repeated with different G operators and the standard deviation between 
the predicted A/3^ values characterize the systematical errors of the matching. Results on larger 
volumes with more blocking levels provide further consistency checks. 



(G(l5;nb,Lb)) = (G((3';nb-l,Lb)) 



(3.1) 



A/3^(/3;7i fo , L b , Oopt) = Ap & (P;n b - l,L h ,a opt ). 



(3.2) 
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Figure 3: Left panel: the bare step scaling function for the Pa/Pf = —0.25 action. The different symbols 
correspond to predictions from optimized matching on different lattice volumes and blocking levels. Right 
panel: the bare step scaling function for the Pa/Pf = —0.15 action. 



Figure 2 illustrates the optimization. The simulations were done on volumes between 32 4 and 
4 4 at P = 6.5 with the Pa/Pf = —0.25 action. The RG blocking was based on 2 HYP smearing 
steps and 5 different operators were used as described in Ref.[8]. 

The left side of the figure shows the optimal A/3 with rib = 3/2/1 blockings and final blocked 
lattices of Lb =2, 3 and 4. The red octagons show the results of the optimized matching, a con- 
sistent value between all three volume series. In contrast the blue diamonds show the predicted 
A/3 without finite volume correction in the optimization. The result on the smallest volume set is 
clearly off, signaling a large finite volume effect. The two larger volumes predict consistent val- 
ues, it appears that with my blocking transformations and 5 operators a finite volume of Lb = 3 is 
already sufficient to minimize finite volume effects. The right side of the figure shows A/3 after 
one more blocking step, with rib = 4/3/2. Again, the finite volume corrected optimized data is sig- 
nificantly different from the uncorrected one but both are consistent with the Lb = 2 results of the 
left hand side. Comparing the finite volume corrected optimized results for A/3 one can conclude 
that A/3 = —0.13(2) on all three volume sequences and after rib = 3/2/1 and ny = 4/3/2 blocking 
levels. The negative value indicates that the RG /3 function has crossed zero, these measurements 
are on the strong coupling side of an IRFP. 

The left panel of Figure 3 shows the result of optimized matching for the Pa/Pf = —0.25 
action at several coupling values, volumes and blocking levels. The data points are consistent 
and show a negative Sb, consistent with the existence of an IRFP at some weaker coupling. MCRG 
matching becomes increasingly sensitive to statistical errors at weaker coupling and I did not pursue 
matching any further with the Pa/Pf = —0.25 action. Rather I considered the action with Pa/Pf = 
—0.15. The result of MCRG matching using rib = 3/2/1 blocking is shown on the right panel of 
Figure 3. This action has a zero in the step scaling function in the investigated region. 
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4. Conclusion 

I have investigated the Nf = 12 flavor SU(3) gauge model. This is a strongly coupled system 
and lattice studies have to control lattice artifacts, in particular the possible appearance of spurious 
fixed points. I showed that the model with plaquette gauge action has such a spurious UVFP in the 
fundamental-adjoint plane that can significantly influence lattice results. Adding a negative adjoint 
plaquette to the gauge action reduces this effect. Using a finite volume improved 2-lattice matching 
MCRG technique I was able to establish the existence of an IRFP for two actions with different 
fundamental-adjoint plaquette ratios, implying infrared conformality for this system. 
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